Halcon查找圆并计算圆的半径

2023-04-06 19:52| 来源: 网络整理| 查看: 265

Halcon查找圆并计算圆的半径原创

茗君（Major_S） 2021-09-07 09:22:29 ©著作权

文章标签 halcon 阈值分割 d3 拟合灰度图文章分类代码人生

©著作权归作者所有：来自51CTO博客作者茗君（Major_S）的原创作品，请联系作者获取转载授权，否则将追究法律责任查找圆并计算圆的半径查找圆并计算圆的半径查找圆并计算圆的半径 connection() select_obj() select_shae() select_shape_xld() segment_contours_xld() select_shape_xld() shape_trans() read_image (Image, 'pumpe') * 1、阈值分割 * threshold (Image, Region, 0, 128) local_threshold (Image, Region, 'adapted_std_deviation', 'dark', [], []) connection (Region, ConnectedRegions) * 2、提取圆形轮廓：select_shape、shape_trans * edges_sub_pix (ConnectedRegions, Edges, 'canny', 1, 40, 70) * select_shape (ConnectedRegions, SelectedRegions, 'circularity', 'and', 1, 100) select_shape (ConnectedRegions, SelectedRegions2, 'outer_radius', 'and', 13, 1000) select_shape (SelectedRegions2, SelectedRegions1, 'circularity', 'and', 0.06, 1) select_shape (SelectedRegions1, SelectedRegions3, 'area', 'and', 75, 99999) shape_trans( SelectedRegions3, RegionTrans, 'outer_circle') dilation_circle (RegionTrans, RegionErosion, 3.5) union1 (RegionErosion, RegionUnion) reduce_domain (Image, RegionUnion, ImageReduced1) connection (ImageReduced1, ConnectedRegions1) count_obj (ConnectedRegions1, Number) Tuple := [] for Index := 1 to Number by 1 select_obj (ConnectedRegions1, ObjectSelected, Index) reduce_domain (ImageReduced1, ObjectSelected, ImageReduced2) threshold_sub_pix (ImageReduced2, Border, 100) select_shape_xld (Border, SelectedXLD, 'contlength', 'and', 35, 9999) segment_contours_xld (SelectedXLD, ContoursSplit, 'lines_circles', 5, 4, 2) fit_circle_contour_xld (SelectedXLD, 'geometric', 100, 0, 0, 8, 2, Row, Column, Radius, StartPhi, EndPhi, PointOrder) Tuple[Index-1]:=Radius[0] endfor stop() * 1.读取图片 read_image (Image, 'pumpe') * 2.创建模板 * 3.匹配模板 * 4.找出区域 * 5.画出轮廓 * 6.计算半径 * 1、阈值分割 * threshold (Image, Region, 0, 128) local_threshold (Image, Region, 'adapted_std_deviation', 'dark', [], []) connection (Region, ConnectedRegions) * 2、提取圆形轮廓：select_shape、shape_trans * edges_sub_pix (ConnectedRegions, Edges, 'canny', 1, 40, 70) * select_shape (ConnectedRegions, SelectedRegions, 'circularity', 'and', 1, 100) select_shape (ConnectedRegions, SelectedRegions, ['outer_radius','anisometry'], 'and',[16,1] ,[30,5]) * 3、从圆形轮廓中提取大圆：dilation_circle dilation_circle (SelectedRegions, RegionDilation, 3.5) * 4、把大圆的灰度图分割出来：reduce_domain union1 (RegionDilation, RegionUnion) reduce_domain (Image, RegionUnion, ImageReduced) edges_sub_pix (ImageReduced, Edges, 'canny', 2, 4, 60) select_shape_xld(Edges, SelectedXLD1, 'area', 'and', 0, 520000) fit_circle_contour_xld(SelectedXLD1, 'algebraic', -1, 0, 0, 3, 2, Row, Column, Radius, StartPhi, EndPhi, PointOrder) gen_circle_contour_xld (ContCircle, Row, Column, Radius, 0, rad(360), 'positive', 1.0) *disp_message(3600, '内轮廓圆心坐标：X=' + Row + ',Y=' + Column +',半径=' + Radius +'.', 'window', -1, -1, 'black', 'true') * 5、从分割出来的图像中提取大圆的轮廓：threshold_sub_pix、select_shape_xld、segment_contours_xld * 6、拟合大圆：fit_circle_contour_xld read_image (Image, 'pumpe') * 1、阈值分割 * threshold (Image, Region, 0, 128) local_threshold (Image, Region, 'adapted_std_deviation', 'dark', [], []) connection (Region, ConnectedRegions) * 2、提取圆形轮廓：select_shape、shape_trans * edges_sub_pix (ConnectedRegions, Edges, 'canny', 1, 40, 70) * select_shape (ConnectedRegions, SelectedRegions, 'circularity', 'and', 1, 100) select_shape (ConnectedRegions, SelectedRegions2, 'outer_radius', 'and', 13, 1000) select_shape (SelectedRegions2, SelectedRegions1, 'circularity', 'and', 0.06, 1) select_shape (SelectedRegions1, SelectedRegions3, 'area', 'and', 75, 99999) shape_trans( SelectedRegions3, RegionTrans, 'outer_circle') dilation_circle (RegionTrans, RegionErosion, 2) union1 (RegionErosion, RegionUnion) reduce_domain (Image, RegionUnion, ImageReduced1) connection (ImageReduced1, ConnectedRegions1) count_obj (ConnectedRegions1, Number) Tuple := [] for Index := 1 to Number by 1 select_obj (ConnectedRegions1, ObjectSelected, Index) reduce_domain (ImageReduced1, ObjectSelected, ImageReduced2) edges_sub_pix(ImageReduced2,Border,'canny',1.1,5,18) *select_shape_xld (Border, SelectedXLD, 'outer_radius', 'and', 17, 22) select_shape_xld (Border, SelectedXLD1, 'contlength', 'and', 30,9000) select_shape_xld (SelectedXLD1, SelectedXLD2, 'outer_radius', 'and', 15.3, 22) fit_circle_contour_xld (SelectedXLD2, 'geometric', 100, 0, 0, 8, 2, Row, Column, Radius, StartPhi, EndPhi, PointOrder) tuple_max (Radius, Max) Tuple[Index-1]:=Max endfor * This program gets the deviation between * circular contours and their best fitting circles * dev_update_var ('off') dev_update_off () * * step: acquire image and initialize the visualization * read_image (Image, 'pumpe') get_image_size (Image, Width, Height) dev_close_window () dev_open_window (0, 0, Width, Height, 'light gray', WindowID) dev_set_part (0, 0, Height - 1, Width - 1) dev_set_line_width (1) dev_set_color ('red') dev_set_draw ('margin') dev_display (Image) set_display_font (WindowID, 16, 'mono', 'true', 'false') disp_continue_message (WindowID, 'black', 'true') stop () * * step: create ROI * fast_threshold (Image, Region, 0, 70, 150) dev_set_colored (3) connection (Region, ConnectedRegions) select_shape (ConnectedRegions, SelectedRegions, ['outer_radius','anisometry','area'], 'and', [5,1,100], [50,1.8,99999]) shape_trans (SelectedRegions, RegionTrans, 'outer_circle') dilation_circle (RegionTrans, RegionDilation, 5.5) union1 (RegionDilation, RegionUnion) reduce_domain (Image, RegionUnion, ImageReduced) dev_clear_window () dev_display (ImageReduced) disp_continue_message (WindowID, 'black', 'true') stop () * * step: create and select relevant contours * threshold_sub_pix (ImageReduced, Border, 80) select_shape_xld (Border, SelectedXLD, ['contlength','outer_radius'], 'and', [70,15], [99999,99999]) segment_contours_xld (SelectedXLD, ContoursSplit, 'lines_circles', 4, 2, 2) select_shape_xld (ContoursSplit, SelectedXLD3, ['outer_radius','contlength'], 'and', [15,30], [45,99999]) union_cocircular_contours_xld (SelectedXLD3, UnionContours2, 0.5, 0.1, 0.2, 2, 10, 10, 'true', 1) sort_contours_xld (UnionContours2, SortedContours, 'upper_left', 'true', 'column') dev_clear_window () dev_set_color ('white') dev_display (ImageReduced) dev_display (SortedContours) disp_continue_message (WindowID, 'black', 'true') stop () * * step: fit circles into the contours and get the * average deviation between contour and circle * count_obj (SortedContours, NumSegments) dev_display (Image) dev_display (SortedContours) NumCircles := 0 disp_message (WindowID, 'Circle radius and average distance', 'window', 10, 10, 'white', 'false') disp_message (WindowID, 'between circle and contour:', 'window', 30, 10, 'white', 'false') for i := 1 to NumSegments by 1 select_obj (SortedContours, SingleSegment, i) NumCircles := NumCircles + 1 fit_circle_contour_xld (SingleSegment, 'atukey', -1, 2, 0, 5, 2, Row, Column, Radius, StartPhi, EndPhi, PointOrder) gen_circle_contour_xld (ContCircle, Row, Column, Radius, 0, rad(360), 'positive', 1) dev_display (ContCircle) dist_ellipse_contour_xld (SingleSegment, 'algebraic', -1, 0, Row, Column, 0, Radius, Radius, MinDist, MaxDist, AvgDist, SigmaDist) disp_message (WindowID, NumCircles, 'window', Row - 10, Column - 5, 'white', 'false') disp_message (WindowID, 'R' + NumCircles + ': ' + Radius$'.3', 'window', (i - 1) * 50 + 30, 450, 'white', 'false') disp_message (WindowID, 'D_avg: ' + AvgDist$'.3', 'window', ((i - 1) * 50) + 50, 450, 'white', 'false') endfor * dev_update_window ('on')

赞收藏评论分享举报

上一篇：halcon定位匹配圆

下一篇：面向对象编程（OOP）六大原则

【本文地址】

公司简介

联系我们

今日新闻

推荐新闻

专题文章